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The dissipative dynamics of a pointlike Josephson junction in binary Bose-condensed mixtures is 
analyzed within the framework of the model of a tunneling Hamiltonian. The transmission of unlike 
particles across a junction is described by the different transmission amplitudes. The effective action 
that describes the dynamics of the phase differences across the junction for each of two condensed 
components is derived employing the functional integration method. In the low-frequency limit the 
dynamics of a Josephson junction can be described by two coupled equations in terms of the potential 
energy and dissipative Rayleigh function using a mechanical analogy. The interplay between mass 
currents of each mixture component appears in the second-order term in the tunneling amplitudes 
due to interspecies hybridizing interaction. The asymmetric case of the binary mixtures with the 
different concentration and order parameters is considered as well. 

PACS numbers: 03.75.Lm, 67.85.Fg, 74.50+r 
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INTRODUCTION 

Recently, experimental study of multicomponent Bosc- 
Einstein condensates has made a substantial progress. 
The study of multiple atomic condensates is intriguing 
since they can produce a laboratory mixture of distin- 
guishable boson superfluids at sufficiently low temper- 
atures. A considerable amount of theoretical work has 
been devoted to binary Bose-condensed mixtures, focus- 
ing, mainly, on the mean-field description of trapped bi- 
nary mixtures [H, stability and phase separation 
collective excitations 0-0] , condensate depletions M, and 
quantum merging of two different condensates [9Mll| . 

In this connection we will consider here the dissipative 
and interference aspects of the Josephson effect in the 
binary Bose-condensed mixtures. The Josephson effect, 
first predicted and discovered for two superconductors 
separated with a thin insulator layer, is a macroscopic 
quantum phenomenon in a condensed medium. The dy- 
namics of the effect is described in terms of the difference 
between the phases of the superconductors, playing the 
role of a macroscopic quantum variable. In spite of wide 
application of the effect in devices for extremely high- 
sensitivity measurements of currents, voltages, and mag- 
netic fields the Josephson effect is still of interest in the 
fundamental modern physics. Like superconductors and 
Fermi superfluids, the Josephson effect is also inherent 
in Bose superfluids [l2l - [l4| . The effect has been observed 
by the mass flow of superfluid 4 He through nanoscale 
apertures coupling two bulk superfluid reservoirs. The 
Josephson tunneling junction in ultracold dilute atomic 
gases is formed with a laser separating two Bose-Einstein 
condensates HI, 

In general, the Josephson effect can include both the 
so-called internal effect for the atoms in different hyper- 
fine states much as NMR phenomena in 3 He pj} and the 
conventional case of two Bose condensates separated with 
a potential barrier which acts as a tunneling junction. A 
lot of work 18-331 has been done in the latter case as a 



direct analogy with conventional superconductors. Those 
studies dealt with the one-component Bose-Einstein con- 
densates alone. 

Multicomponent Bose-Einstein condensates are also a 
very interesting subject for studying various macroscopic 
tunneling phenomena (33Tl37j . One may expect novel and 
richer manifestations of the Josephson effect. For the 
system of two Bose-condensed mixtures connected with 
a weakly coupled junction, the dynamics of the Joseph- 
son effect should be governed by the difference between 
the phases for each Bose-condensed component of a mix- 
ture. In other words, two relative phases ipi and cp2 must 
be involved into consideration. In addition, we must take 
into account the different tunneling transition amplitudes 
Ji and I2 across the junction for various bosonic atoms 
of masses mi and m-i composing the mixture. This re- 
sults in two Josephson currents associated with the mass 
flow of each component of a mixture across the junction. 
From the general point of view one may expect interfer- 
ence and coupling between the Josephson currents. The 
dynamics of the Josephson effect in the vicinity of the 
phase separation of a mixture is an additional motiva- 
tion for studying binary condensates. Furthermore, the 
dissipative aspects of the Josephson dynamics in binary 
condensed mixtures have not yet received a proper and 
wide investigation. 

The dissipative effects and dephasing of the Josephson 
oscillations come from the coupling between the macro- 
scopic relative phase variable and the infinite number 
of the microscopic degrees of freedom. The successive 
method of eliminating microscopic degrees of freedom 
from the Hamiltonian was developed first for the su- 
perconducting Josephson systems Sill. Later that 
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functional integration approach 
Green function method [301 ] were extended and applied 
to studying dissipative and nonequilibrium Josephson dy- 
namics in the one-component Bose-condensed systems. 

In this paper, we will generalize the energy dissipa- 
tion effects in the Josephson dynamics to the case of the 
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binary Bose-condensed systems, employing functional in- 
tegration approach of Ref. [28]. We will derive the ex- 
pression for the effective action depending on two relative 
phases ifi and (p 2 between two condensed mixtures con- 
nected by a pointlike tunneling junction. The response 
functions in the effective action give the full information 
on the dynamics of the junction. The low-frequency ex- 
pansion of the response functions allows us to determine 
two coupled Josephson equations for the relative phases 
f i and tp 2 , Josephson energy U(cp\, ip 2 ) and dissipative 
Rayleigh function R(<px, <j>2)- Of course, we consider the 
region of the parameters in which the homogeneous state 
of the both mixtures in the left-hand and right-hand 
bulks is stable and the mixtures are not phase-separated. 



EFFECTIVE ACTION 

We keep in mind the case of a pointlike and weakly 
coupled junction between two macroscopic infinite reser- 
voirs containing binary condensed mixtures. In addition, 
we neglect the feedback effect of the junction on the mix- 
tures and assume that both the mixtures are always in 
the thermal equilibrium state. The image of the sys- 
tem is two bulks with one common point through which 
the transmission of particles is possible with the different 
tunneling amplitudes I\ and I 2 depending on the type of 
particles. 

So, our starting point is the so-called tunneling Hamil- 
tonian (h = 1, volume V = 1) 



H = H l +H r + H u + H t , 



(1) 



where Hi r describes the bulk binary Bose-condensed 
mixture on the left-hand and right-hand sides, respec- 
tively, 

+\ ]T u ik ,i [d 3 r^l^l ,,¥ M . 



i,fc=l,2 

Here i and k take 1 or 2 and denote one of the com- 
ponents of a mixture composed with particles of mass 
mi and m 2 . The coupling between particles is specified 
by the constants Uik, 1 = Uki, 1 which can be expressed 
by means of the s-scattering length a^, ; according to 
Uik,i — 2™^, ((m^ 1 + m~ k 1 ). The same expressions with 
the substitution I — > r refer to the mixture on the right- 
hand side. 

The energy, associated with varying the number of par- 
ticles on the left-hand and right-hand sides, 
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Hv = - 



i, k=l, 2 



N, 



i.l 



N, 



i,r TT Nk,l 
Uik 



is analogous to the capacity energy of a junction in the 
case of superconductors. The constants Uik can be con- 
nected with the second derivatives of the total energy 

E = E[(N 1J , N 2J ), (JVi, r , N 2 , r )] 
= E l {N hh N 2 ,i) + E r {N hr> N 2 , r ) 

with respect to the relative change in the number of the 
particles across the junction, 



Uik 



d 2 E 



+ 



d 2 E 



dNijdNkj dN^ r dN k , r 
d 2 E l d 2 E r 



dNi,idN k ,i dNi^dNk^' 

with the obvious symmetrical relation Uik — Uki- 
constants Uik can usually be estimated also as 



The 



U, 



ik 



dfMi 



dN k 



k, 1 



dN, 



k, r 



The term Hjj describes the point that the energy of the 
system on the whole may depend on the relative numbers 
of particles from the left-hand and right-hand bulks. In 
order to avoid an instability of the total system against 
an infinite growth of the number of particles on the left- 
hand or right-hand sides, it is necessary to suppose that 
Hij > for any variations Ni t 1 — Ni t r and N% / — N 2 . r - 
In other words, the matrix of coefficients Uik should be 
positively determined, i.e., ||f7ft;|| > 0. The total number 
of the particles of the type labelled by i in each bulk I or 
r is given by 



Nij= ld a r^l^ij 



and (I — > r) . 



The last term, 
H t = -f d 3 rd 3 r' [Vl ^hir, r')*i,r(r') 

+*5 )I (r)7 2 (r, r')*2,r(r')+H.c], 

is responsible for the transitions of particles from the 
right-hand to the left-hand bulk and vice versa. In gen- 
eral, the transition amplitudes I± and I 2 are different for 
the various species of the particles composing the mix- 
ture. We consider here the simplest case of a pointlikc 
junction, i.e., 

Ii(r, r') = IiS{r)S{r'), i = l,2. 

To study the properties of the system described by 
Eq. we calculate the partition function Z using the 
analogy of the superconducting junction 3|*; 4 0| an d the 
approaches employed for the Bose junction [271. |28|. 

Z = J 2? 2 *i,,2? 2 *2,l2? 2 *l,r2? 2 *2,rex P [-5B], 
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where the action Se, defined on the imaginary (Matsub- 
ara) time r, reads 

S E = cLtLe, 

J-P/2 

i, k=l, 2 X / 

and, as usual, j3 = 1/T is an inverse temperature. 

To eliminate the quartic terms in the action, which 
come from the energy Hjj, we employ the Hubbard- 
Stratonovich procedure by introducing additional gauge 
fields Vi (t) and V 2 (r) in analogy with the so-called plas- 
mon gauge field in metals, 

exp(-J drHu^j = Jw 1 (t)VV 2 (t) 



exp 



N,,-N, 



z,k=l,2 



J VVi (r) VV 2 (r) cxp - J dr ^ 
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1. 



i,fc=l,2 



Here U^ 1 is an inverse matrix for the 2x2 matrix Uik 



r j-i = 1 / U22 -u 12 

lh UnU 22 - U 12 U 2 i \ -U21 U u 



After the introduction of the fields V\ (r) and V 2 (r) the 
partition function Z takes the form 

Z= fw 1 VV 2 [] P 2 * i;i P 2 * j;r exp[-5 B ], 

^ »=1,2 

where 

S E = Si + S r + I dTH t + \ J2 I drViU^Vk. 

J i.k=1.2 J 



Here 5; and S r denote 



i=l,2 



^ = 5 --^E dTd 3 r^l r V t ( T )^ hr . 



i=l,2 

In essence, this replacement looks like the renormaliza- 
tion of the chemical potentials for each component of a 
mixture on the left and right-hand sides of the junction, 

MM -> IH,l - iVi{r)/2 , 
Mi,r -> A*i,r + i^i(r)/2 . 

At this point it is advantageous to perform a gauge 
transformation of the field operators ^ and "J, which 



makes the future Green functions real. This is achieved 
by introducing the phases ip^ 1 and ipi t r according to 

*», 1 (», r) -> exp [i<p it 1 (i) r) (r)] j (ii r) , 

exp[-*Vi,I(i.r)( r )] *i,J(i,r)> 



t 

i, I (i, r) 



and by imposing the conditions (p^.; = — Vi(r)/2 and 
<Pj )7 . = Vi(r)/2. Thus, we arrive at the first Josephson 
relations for the relative phases <Pi(t) and <P 2 (t): 

WW = ^iW and <Pi = <Pi, r ~ <Pi,i (i = 1, 2). (2) 

That we have achieved is only a formal elimination of the 
explicit dependence of the chemical potentials /x^ ; and 
r upon the time r. On the other hand, the tunneling 
amplitudes 7i and I 2 acquire additional factors depend- 
ing on the phase differences (pi(r) and ip 2 (r) across the 
junction, 

7i -> A = Jie'^M and J 2 -> J 2 = hc lip ^ T) . 
Hence we arrive at 



=1,2 



V = <,'„ + fdrH t + I dr 

J J „- 7 1 



i, fe=l, 2 

where S = 5; + SV- The tunneling term is given by 
Ht = - j e^rd 3 / 2 [*t J J i (p,p'; T )* < , I . + H.c.]. 



Next, one must integrate over fields ^ and W in order 
to obtain the effective action S e s depending on V\ (r) and 
V 2 (t) alone. In calculations we treat H t as a perturba- 
tion and restrict ourselves by second-order perturbation 
in the tunneling amplitudes I\ and I 2 . Omitting the term 
independent of Vi(r) and employing the Josephson rela- 
tions Vj(r) = <pi(r), we find the effective action as 



((H?))o 



S e s[(Pi,<P2] = I dry Y Vi-^^k + (H t ) - 

i,fe=l,2 



Here (A)q means the averaging over decoupled action 
So = Si + S r corresponding to Hq = Hi + H r , i.e., 

= ( ^!^ 0) and «ff 2 » = (Hf) a (H t f . 

It is obvious that 

(Ht)o = (Hit)o + (H 2 t)o 
((Hf)) = ((Hf t )) + 2((H lt H 2t )) + «fff t )) . 

The first-order terms in the tunneling transparency are 
obviously decoupled. Since Hi t H 2t <~ /1/2, second-order 
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terms will result in the coupling and interference between 
the mass currents of atom species 1 and 2. Averaging over 
the I and r variables is independent of each other. 

In the course of calculation we follow the Bogoliubov 
method of separating the field operators into the conden- 
sate G and noncondensate $ fractions, i.e., 

^i, /(z, r) t) t) > 

with the conventional relation Cj 5 r ) — n it r j 



where ni t iu r \ is the density of particles labelled with 
i = 1, 2 in the condensate fraction in the left-hand and 
right-hand bulks, respectively. 

In a binary Bose-condensed mixture the Green func- 
tion G(uj n ,p) represents a block 4x4 matrix. The Green 
function can readily be found from the inverse matrix 
whose Fourier representation in the approximation of a 
weakly interacting two-component Bose-condensed gas 
mixture is given by 



G- l (u> n ,p)= 



_/G n 1 A 12 



A12 G 2 



+ m 

A n 

A12 
A12 



A n 



An 

- i]i + An 

A12 
A12 



Ai 2 
A12 
+ V2 

A 22 



\ 22 



A12 
A12 

A 22 

h T] 2 + A 2 2/ 



Here w„ = lixnT is the Matsubara frequency and p is 
the momentum. Also we have introduced the following 
notations for the free-particle energies 

m = Vi(p) = p 2 /2m 1 , ?]2 = 772 (p) = p 2 /2m 2 ; 

and for the order parameters 

An=uiini, A 22 = U22W2 , A12 = u 12 ^nin 2 . 



Accordingly, for the direct matrix Green function 



G(u n , p) = 



G\i G12 
G21 G22 



I Gn F\\ G12 F\2 \ 

F\x Gn F\ 2 Gyi 

G21 Fix G22 F22 

\F\ X G21 F\ 2 G22) 



we arrive at the following components of the matrix: 



Gn(cj„, p) = G n (-u„, p) = 



(iuj n + 771 + A u )(^ + e 2 2 ) - 2A 2 2 ?7 2 



Fu(u n , p) = F n (-w„, p) 
Gi 2 (w„, p) = G\2{-u) n , p) 
F12 p) = F\ 2 {~uj ni p) = 



-An(^ + £ i) + 2A^ 2 

(w 2 +W 2 )( W 2 +W 2 ) 
_ Ai2^„ + T]i){iuj n + 772) 

K + w 2 )K + c 2 2 ) 

A 12 {iuj n + r]i)(-ioj n + r? 2 ) 

( W 2 + W 2 )(^ 2 + W 2 ) 



The lower part of the Green function matrix is deter- 
mined by permutation 1^2 and A21 = A12. We have 
introduced above the similar abbreviations for either of 
two mixtures, using the Bogoliubov nomenclature 
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vi 



2An?7i and e 2 



2A 22 ??2 



where 77^ 2 = p 1 j2m\ 2 is the free-particle energy. The 
energies S\ and £2 are the Bogoliubov energies of each 
component of a mixture taken separately. In the mix- 
ture the interspecies interaction u\2 hybridizes these two 
modes, resulting in two familiar branches of elementary 



excitation spectrum, e.g., 3, 8] 



;?, 2 (p) = §(ei 



4± 



(ej - e 2 ) 2 + l6Al 2mm 



with the crossover to the sound-like dispersion wi. 2 = 
pc\,2 at small p — > momentum. The sound velocities 
ci 2 in a mixture are determined by the well-known rela- 
tions as well d, Q 



c l,2 



lrAn 



All + A 2 2 ± //An A 22 
777,1 771 2 \/ V 777,1 777,2 



Ai 

77*1 



A 

m 2 



A 2 
mi?n 2 
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The inequalities u\\u 22 > u\ 2 and c\ 2 > are certainly 
supposed to guarantee the stability of a mixture against 
its dcmixing. 

Taking into account @ , we arrive finally at the follow- 
ing generalization of the effective action compared with 
that in the one-component condensed system [28} 



S e g[ipi(T),cp 2 (T)} = / dr [ ^ <&( 

^ i,fc=l,2 



— 2Ily/nuni r COS^l(r) - 2I 2y /n 2 in 2r COS(/? 2 (t) 



i,fc=l,2 



- r') cos ^i(r) - <Pfc(r')) 



+Afc(r - t') cosher) + (pk(r') 



(3) 



Here and /3i k are the so-called response functions 
which can be written using the Green functions: 

atk(r) = V n n n kigtk,r( T ) + \ /n irn k rgtkj( T ) + Gik(r), 
Pik{r) = y/nuriki ftk, r ( T ) + yjn-irrikr ft k ,i( T ) + J 7 ik{r). 

The Green functions are calculated at the junction point, 
i.e., at r = and r' = 0, 

9tk,l(r)( T ) = \ \9ikJ(r){T) + 9ik,l{r){r)\ 

d 3 p 



2(2tt) 



[G ik , i( r ) (p, t) + G lfei l(r ) (p, r)] , 



fik,l(r)( T ) 



d 3 p 



[F ik<l{r) (p, t) + F*, r) (p, r)]: 



2(27r) 3 rws.HT-jv-f •) 1 i 4 fe,i(r)^ 

0ifc(V) = 9ik,l{T)gik,r{T) + Sik,l(T)gik,r(T), 
Fik(T) = /ift,i( T )/iife,r( r ) + fL,l( T )f^,r(r). 

In order to comprehend the dynamics of the rela- 
tive phase differences tpi and tp 2 across the junction, we 
should analyze the behavior of the response functions 
ctikir) and /3jfc(r) as a function of time. Note that the 
contribution of the terms Gik and Tik to the response 
functions a%k and /3jk is much smaller than that of the 
first two others 28j. The order-of- magnitude smallness 
is about a ratio of the noncondensate density to the con- 
densate density or about gas parameter (na 3 ) 1 / 2 <C 1. 
Below, analyzing and ft^, we will concentrate our 
attention on the first two terms which can be attributed 
to the condensate-noncondensate tunneling processes. 



THE RESPONSE FUNCTIONS 

The calculation of the response functions in the gen- 
eral form in a mixture is a complicated problem. Keep- 
ing in mind the study of the low-frequency dynamics of a 
junction, we will restrict our calculation by analyzing the 
behavior of the response functions on the long-time scale. 



This means that we should find the low-frequency decom- 
position of the response functions in the Matsubara fre- 
quencies u n . In fact, we imply the inequality \uj n \ <C u)\, 
lu 2 . Next, we will employ the procedure of analytical 
continuation in order to derive the dynamic Josephson 
equations which the relative phases <p\ and tp 2 obey. To 
obtain the dissipative terms, it is sufficient to expand 
the response functions (Xik and fiik up to linear terms in 
frequency uj„. 

So, we look for the following first coefficients in the 
low-frequency decomposition 



* (0) 



Mi 



aik{u n ) = —c 

Pik(Un) = -P& ) +P$\u>n\+-~ 

Accordingly, the expressions for the response functions 
in the imaginary-time representation read as 

2 



fe(r) = -4° ) ^)-4 1) ^G 



1)1/ *T V 
fe 7T Vsin(7rTT)/ 



+ ■ 



7TT 



-V 



vSin(7rTr)/ 

First of all, we should note that zero harmonic u n = 
in the response is unimportant if i = k, and thus 
we can deal with the difference au(ui n ) = au(uj n ) — 
cnii(ui n = 0). In fact, this corresponds to the substitution 
a u( T ) = &u(t) +o;,i(0)5(r) into effective action The 
second term au(0)5(T) yields a physically unimportant 
time- and phase-independent contribution to the action, 
meaning a shift of the ground-state energy of a junction. 
For i ^ k, this does not hold for. As we will see below, 
a i^k an d flik are connected with the Josephson potential 

energy and , determine the dissipative properties 
of the junction. 

Let us start from For the sake of brevity, we 

present here the expressions for aik and /3jfc in the case 
of a symmetric junction with the identical mixtures on 
the both left-hand and right-hand sides of the junction 
when A lk j = A ik ^ r = A lk and mi = n lr = m. The 
general case I 7^ r will be given in the Appendix. The 
simple calculation yields 



-2m/+K = 0) 

m\ (An + m 2 c 1 c 2 ) 



7r (miAn + 2mim 2 cic 2 + m 2 A 22 ) 1 / 2 ' 

$ = -2^/^fU" = 0) = 
mi?7i 2 Ai 2 



P12 



\/nin 2 



7r (miAn + 2mim 2 cic 2 + m 2 A 22 ) 1 / 2 ' 
and the other quantities can be obtained with 1^2 



(4) 



The calculation of and a^} is more complicated: 



„M _ 0(1) _ 



mi cic 2 + A 22 /m 2 



47TCiC 2 

a W = /3 W = 2^h~ 2 



ci + c 2 
1 A12 



47TCiC 2 Ci + C 2 



(5) 
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Again all the remaining quantities are given by 1 2. 
Note only that all fi^) and ocf^ k remain finite and nonsin- 
gular at the demixing point A11A22 = A^ 2 or when one 
of the sound velocities vanishes c 2 = 0. On the contrary 
both and $ 
point as c 2 — » 0. 



both a^jP and diverge with approaching at demixing 



JOSEPHSON EQUATIONS 

To obtain the dynamics of the relative phases tpi and 
ip 2 in real time, we now follow the standard procedure of 
analytical continuation. Accordingly, the substitution of 
Matsubara frequencies \u> n \ — > —iu> in the Fourier trans- 
form of the Euler-Lagrange equations dS e f f /dipi(r) = 0, 

- U^ifkir) + 2Iim sin ^i(r) 

k=l,2 

+27, ^2 !k I dr' aik(r - t') sin(ipi(r) - ipk(T')\ 

1 1 o « 



fc=l,2 



+^ife( / r-T / )sin^ i (r) +^(t')J 



= 0, 



will yield the real-time equation for the phase <fi dynam- 
ics. 

As we have found in the limit of the slowly varying 
phases </?i,2, the response functions in the real-time rep- 
resentations go over to the following decomposition 

a ik (t) = -a^S(t)+a^5'(t) + --- 

Then, we can derive a couple of the Josephson equations 
which every relative phase ipiit) obeys in the real time t: 

^ U^}<pk(t) + 2Iini sin ifi(t) + ^ 27;/ fe x 



fc=l,2 



fe=l,2 



[a$ sin(^(i) - ^ fe (t))+/4 0) Bin( W (t) + p fc (t) 
+ 2IiI k (p k (t)x 



fc=l,2 



^ cos^(t) - + cos^i(t) + 9?fc(t) 



= 0. 



Using relations from (|4]) and ((5]), we finally arrive at the 
desired Josephson equations for two phases <pi(t) and 

</*(!)•• 

^ ^Vfe + 4/i/ fc a^Vfc c os<y5iCOS^ fc 

fe=l,2 fe=l,2 

+21^ sin + ^ 4IiI k /3^ sin ^ cos <p k = 0, i = 1, 2. 



fe=l,2 



We can interpret these Josephson equations using a 
mechanical analogy. First, we introduce the potential 



energy U(<pi, (p 2 ) of a junction according to 
U (<px, ip 2 ) = ~y] cos + i ^ e ifc cos ^ cos 

i=l,2 t,fc=l,2 

with the energy coefficients determined by 

Ei = 2hn t and e« = 47, J^ 5 , (f, fc = 1, 2). 

Then we can introduce the Lagrangian L = K — U as a 
difference between the kinetic K and potential U energies 



L = K-U= J2 



u^ 1 



tptifik - U(<pi, ip 2 )- 



i,fc=l,2 



The dissipation energy effect can be described by intro- 
ducing the dissipative Rayleigh function according to 



i,A=l,2 



(i) 



,fe=l,2 



Finally, the Josephson equations can be written in the 
general form as 

d dL dL dR 
— i = 12, 

dt dipi d<pi d<pi ' 

where in accordance with Eq. ^ , 

(pi(t) = -5fii(t) = Hi, l- lM,r- 

The dissipative function R has a sense of the energy 
dissipation power in the system. This is obvious from the 
following equation, 

dH d /^—v . dL \ x - . dR nr> 



which means that the energy dissipation power equals the 
double dissipative function. Since the energy dissipation 
must result in decreasing the total energy 77 = K + U 
of the system, the dissipative function R must be a pos- 
itively determined matrix R > 0, i.e., 

rn > and rur 22 > r 12 r 2 Y. 

We are persuaded that this is true from 



nir 22 
r\ 2 r 2 \ 



- 1 = 



1 = 



a , a 



12 "21 



mim 2 cic 2 (ci + c 2 ) 5 

^12 



> 0. 



Note that the condition R > gets broken simultane- 
ously with the condition 01,2 > necessary for the sta- 
bility of a mixture against its demixing. In addition, we 
also disclose a symmetry of kinetic dissipative coefficients 
Tik = Tki in accordance with the Onsager principle. 
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It is interesting that the point of demixing instability 
c\C2 = is not singular for the potential energy coeffi- 
cients Eik- On the contrary, the dissipative coefficients 
rjfe become infinite. The latter means that the Joseph- 
son dynamics should slow down and demonstrate an en- 
hancement of decoherence and damping of the Josephson 
oscillations in the vicinity of the phase demixing. From 
the mechanical point of view the low-frequency dynam- 
ics of a Josephson junction in a Bose-condensed mixture 
can be described as a system of two coupled particles or 
pendula moving or oscillating in a viscous medium in a 
periodic potential relief. 

In the lack of hybridization between the different atom 
species, i.e., when A i2 = 0, the crossed terms in the re- 
sponse functions a^k and ft^fc vanish as well. As both 
Oii^k = and Pi^k = 0, the Josephson equations split 
into two decoupled equations for each component of a 
mixture. In this case no interference in the mass cur- 
rents appears and the response functions together with 
nonzero diagonal coefficients an , fin go over to the quan- 
tities corresponding to the case of a single-component 
Bose-condensed gas 27, 2£jj. 



are weakly damped. Emphasize that the interference be- 
tween the Josephson and dissipative Ohmic components 
of a mass current for each type of particles starts only 
from the second-order terms in the tunneling amplitudes 
and, eventually, due to the presence of the noncondensate 
fractions. The interference entails, in particular, that 
the maximum amplitude of the Josephson current of one 
species atoms depends on the relative phase difference of 
the second component of a mixture. In addition, it be- 
comes possible that the imbalance in the chemical poten- 
tial of one component of a mixture can induce also the 
Ohmic contribution into the mass current of the other 
component. We believe these aspects deserve a further 
study. 
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APPENDIX 



CONCLUSION 

To summarize, in this paper we have used a functional 
integration method for the model of a tunneling Hamil- 
tonian in order to analyze the energy dissipation effects 
in the dynamics of a pointlike Josephson junction be- 
tween two weakly nonideal Bose-condensed gas mixtures 
in the thermal equilibrium. The transmission of parti- 
cles of each component of a mixture across the junction 
is described by two different tunneling amplitudes. The 
effective action and response functions that describe the 
dynamics of two relative phases ipi and if2 correspond- 
ing to each condensed component of a mixture are found. 
The quasiclassical Josephson equations for the relative 
phases are derived from the low-frequency decomposition 
of the response functions. 

The dynamics of a pointlike junction in a mixture dis- 
plays a dissipative Ohmic nature. The energy dissipation 
effects result from the noncondensate excitations and ap- 
pear in second order in the tunneling amplitudes. The 
latter fact favors low damping rates of the Josephson os- 
cillations in the pointlike junctions. A growth of the tem- 
perature leads to decreasing the Josephson energy and 
to increasing the energy dissipation power. The close- 
ness to the phase separation of a mixture enhances the 
Ohmic character of the phase dynamics. The dissipative 
Rayleigh function is determined. 

On the whole, the Josephson phase dynamics in bi- 
nary mixtures is described by two coupled equations. 
This means, in particular, an existence of two Josephson 
frequencies for small oscillations of the phases around 
ipi = and ip 2 — 0. Since ^ 0, the oscillations 



Here we present the general asymmetric case of a junc- 
tion if nn ^ riir and A^,/ ^ A^, r . For the Josephson 
potential energy 

U((fil , If 2 ) = - ^ Ei COS Ifi + - ^ £ ik c OS ft COS Ifik, 

we find 



i=l, 2 



. fc=l. 2 



Ei = 2Iiy/minir (i = 1) 2) 
and the next terms 

2lfml r nu(A llir + m 2 ci r c 2r ) 



£ll 



(miA UjT . + 2mim 2 ci r c 2r + m 2 A 2 2,r) 1 / 2 
ni r (An t i + 1TI2CUC21) 



£22 



(miAii.j + 2miTO 2 ci/c 2 ; + m 2 A 22 ,i) 1 / 2 
n-2/(A 22 , r + mici r c 2r ) 



27|m| 



(miAn ir + 2m 1 m,2C lr C2 r + W2A 22 , r ) 1 / 2 
n 2r -(A 22 j + mic u c2i) 



(roiAii,) + 2m 1 m 2 c 1 iC2i + m 2 A 22 ,i) 1/2 



The other two nondiagonal coefficients e\2 = £21 are 
given by the expression 



£12 = 



2/1/2 



V /n l'"2/ Ai 2i r y/n lr n2r A 12 ,/ 



C2r 



CUC21 Cu + C 2 l 



For the kinetic coefficients in the dissipative Rayleigh 
function 



R(fl, <f2) 



\ J2 r ^ 

i=l, 2 



COS Ifi COS (fkflf2- 



8 



wc have 



Tl2 



If mi r n u c lr c 2r + A 2 2, r /m 2 



7T 



Ci r C 2r 



+ 



nir ci/c 2 ; 



c lr + c 2r 

A 22)Z /m 2 " 



CUC21 
2 

2 

7T 



7|m 2 



Cli + C 2 ; 

«2i c lr c 2r + An if ./mi 



c lr c 2r 



Clr T C 2r 



+ 



«2r ci;c 2 ; + An, i /mi 



C11C21 



en + c 2 i 



The other two nondiagonal coefficients n 2 
found from the expression 



r 2 i can be 



ri2 



nun 2 l Ai2, r , ni r n 2r A 



ci r c 2r ci 7 



C2r 



+ 



12J 



Cl/C 2 ; Cil + C 2 l 



The expressions derived above for the dynamical coef- 
ficients in the Josephson equations governing the phase 
differences (pi(t) and (p 2 (t) across a pointlike junction 
allow us to describe the low-frequency dynamics in the 
asymmetric case of Bose-condensed gas mixtures with the 
different densities of the atom species and with the differ- 
ent order parameters. The stability of mixtures against 
phase separation and the positive definiteness of the dis- 
sipative Rayleigh function imply the fulfillment of condi- 
tions en, C21 > and ci r , c 2r > 0. 
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